
//-------------------------------
// FIGURE 1
//-------------------------------

clear
use masterdata_policing.dta


local Y lfinesandforfeitspcp

forvalues i = 1(1)1000 {
	local log`i' = log(`i'+1)
}
local loghalf = log(.5 + 1)

*** do the graph ***

su `Y', det

#delimit;

hist `Y'
	,
	xlab(
		0 "0"
		`loghalf' ".5"
		`log1' "1"
		`log2' "2"
		`log5' "5"
		`log10' "10"
		`log20' "20"
		`log35' "35"
		`log100' "100"
		`log500' "500"
		,
			angle(horiz)
			
	)
	xline(`r(mean)', lcol(gray) lpat(dash) lwid(thick))
	fcol(none)
	//lwid(medium)
	lcol(black)
	//barwidth(.08)
	graphregion(color(white)) 
	plotregion(style(none))
	xtitle("Fine revenue per capita")
	ytitle("Number of cities")
	freq
	ylab(, angle(horiz))
	saving(hist2, replace)
	;
		
#delimit cr

gr combine hist1.gph hist2.gph, ysize(3) graphregion(color(white)) 

